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Abstract. When an electromagnetic signal propagates in vacuo, a polarization detec- 
tor cannot be rigorously perpendicular to the wave vector because of diffraction effects. 
The vacuum behaves as a noisy channel, even if the detectors are perfect. The "noise" 
can however be reduced and nearly cancelled by a relative motion of the observer toward 
the source. The standard definition of a reduced density matrix fails for photon polar- 
ization, because the transversality condition behaves like a superselection rule. We can 
however define an effective reduced density matrix which corresponds to a restricted class 
of positive operator-valued measures. There are no pure photon qubits, and no exactly 
orthogonal qubit states. 



1. Introduction 

The long range propagation of polarized photons is an essential tool of quantum crypto- 
graphy [1]. Usually, optical fibers are used, and the photons may be absorbed or depo- 
larized due to imperfections. In some cases, such as communication with space stations, 
the photons must propagate in vacuo [2]. The beam then has a finite diffraction angle of 
order A/a, where a is the aperture size, and new deleterious effects appear. In particular 
a polarization detector cannot be rigorously perpendicular to the wave vector, and the 
transmission is never faithful, even with perfect detectors. Moreover, the "vacuum noise" 
depends on the relative motion of the observer with respect to the source. 

The relativistic effects reported here are essentially different from those for massive 
particles [3] because massless particles have only two linearly independent polarization 
states. The properties that we discuss are kinematical, not dynamical. At the statistical 
level, it is not even necessary to involve quantum electrodynamics. Most formulas can 
be derived by elementary classical methods as shown below. It is only when we need to 
consider individual photons, for cryptographic applications, that quantum theory becomes 
essential. 

This article consists of two parts. First we consider the propagation of a classical 
electromagnetic wave. The wave vector cannot be constant because of diffraction effects. 
A polarization detector cannot unambiguously distinguish orthogonal polarizations, even 
if the detector is perfect. The vacuum behaves as a noisy channel. We then show that this 
"noise" can be reduced and nearly cancelled by a relative motion of the observer toward 
the source. 

In the second part of this paper, we investigate the transmission of a single photon. The 
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diffraction effects mentioned above lead to superselection rules which make it impossible 
to define a reduced density matrix for polarization. It is still possible to have "effective" 
density matrices; however, the latter depend not only on the preparation process, but also 
on the type of detection that is used by the observer. 

2. Classical electromagnetic signals 

Assume for simplicity that the electromagnetic signal is monochromatic. In a Fourier 
decomposition, the Cartesian components of the wave vector (with /i — 0, 1,2,3) can 
be written in term of polar angles: 

/c M = (1, sin 9 cos 0, sin 9 sin 0, cos 9), (1) 

where we use units such that c = 1 and k = 1. Let us choose the z axis so that a well 
collimated beam has a large amplitude only for small 9, and let us rotate the x and y 
axes so that a particular k in which we are interested has = (we shall later return to 
arbitrary k with 0^0). The Fourier transform of the electric field is perpendicular to 



If the emitter ("Alice") selects a linear polarization angle a, the Fourier transform of E 
is proportional to 



The magnetic field Fourier transform is proportional to B k = k x E k . The Poynting 
vector P is parallel to k and gives the energy flux, as usual. We shall henceforth omit the 
subscript k. 

Suppose that the receiver ("Bob") has an infinite flat detector parallel to the xy plane. 
Then only the component 

P z oc E 2 X + Ey = cos 2 a cos 2 9 + sin 2 a (4) 

is absorbed by the detector. Moreover, if Bob selects a polarization angle b in the xy 
plane, the flux detected by him will be proportional to 

(E x cos b + E y sin b) 2 = (cos a cos b cos 9 + sin a sin b) 2 . (5) 

Under ideal conditions (9 = 0), the maximal signal behaves as cos 2 (a — b). For small but 
finite 9, the fraction of the signal that is lost, when a = b, is 9 2 cos 2 a. More generally it is 
9 2 cos 2 (a — 0), where is the azimuthal angle of k, which was set to zero in the preceding 
calculation by a suitable choice of the x axis. 

Apart from the above loss in intensity, the angle mismatch may introduce errors. Con- 
sider again the case = 0. When Alice emits a wave that is linearly polarized along 
the x or y axes, Bob's intensity pattern varies as A 2 cos 2 b and A 2 sin 2 b, respectively. 
These intensity distributions would be unambiguously recognized as pertaining to the 
corresponding linear polarizations. However, for a general linear polarization state that 
Alice may send (0 ^ a ^ tt/2) the intensity will be distributed according to Eq. (5). As a 
result, the angle a' that Bob will ascribe to it would be related to the actual polarization 
direction by tana' = tana/ cos 9. 



k = (sin#, 0, cos#). 



(2) 



Ek = (cos a cos 9, sin a, — cos a sin 9) . 



(3) 
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In a real experiment, the angles 6 and are distributed in a continuous way around 
the z axis (exactly how depends on the properties of the laser) and one has to take a 
suitable average over them. Since the definition of polarization explicitly depends on the 
direction of k, taking the average over many values of k leads to an impure polarization 
and therefore may cause not only an attenuation of the beam, but also identification 
errors. 

Let us now consider the effect of a motion of Bob relative to Alice, with a constant 
velocity v = (0, 0, v). The Lorentz transformation of in Eq. (1) yields new components 

k' = 7(1 — v cos 6) and k' z = 7(cos6> — v), (6) 

where 7 = (1 — v 2 )~ 1 / 2 . Considering again a single Fourier component, we have, instead 
of the unit vector k, a new unit vector 

k , = / sine ^ cosO-v \ 
\7(1 — v cos6>) ' ' l — vcos6J 

In other words, there is a new tilt angle 6' given by 

sin#' = sin 0/7(1 - vco$6). (8) 

For small 8, such that 6 2 <C \v |, we have 




d ' = d \h—, (9) 



The square root is the familiar relativistic Doppler factor. For large negative v, the 
diffraction angle becomes arbitrarily small, and sideway losses can be reduced to zero. 

It is noteworthy that the same Doppler factor was obtained by Jarett and Cover [4] who 
considered only the relativistic transformations of bit rate and noise intensity, without 
any specific physical model. This remarkable agreement shows that information theory 
should properly be considered as a branch of physics. 

3. Quantized electromagnetic signals 

We now turn to the quantum language in order to discuss applications to secure com- 
munication. The ideal scenario is that Alice sends isolated photons (one particle Fock 
states). In a more realistic setup, the transmission is by means of weak coherent pulses 
containing on the average less than one photon each. 

A basis of the one-photon space is spanned by states of definite momentum and helicity, 

|k,e±) = |k)®|e±>, (10) 

where helicity states |e^) are explicitly defined in Eq. (15) below. The momentum basis 
is normalized by 

(q|k) = (27r) 3 (2A;V 3) (q-k). (11) 

Polarization states that correspond to different momenta belong to distinct Hilbert 
spaces and cannot be superposed (an expression such as |e^) + |€q) is meaningless if 
k 7^ q). The complete basis (10) does not violate this superselection rule, owing to the 
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othogonality of the momentum basis. Therefore, a generic one-photon state is given by a 
wave packet [5] 

|tf) = /^(k)/(k)|k,a(k)>, (12) 

where the polarization state |a(k)) corresponds to the 3- vector 

a(k) = a + (k)e+ + «_(k)e k , (13) 

|a+| 2 + |«-| 2 = 1, and the explicit form of e k is given below. The Lorentz-invariant 
measure is 

* w = (Sow' (14) 

and normalized states satisfy / dfi(k)\f(\t)\ 2 = 1. Since diffraction angles and frequency 
spreads are usually small [2, 5], /(k) significantly differs from zero only in the vicinity of 
a certain momentum that we shall denote by k^. 

Lorentz transformations of quantum states are most easily computed by referring to 
some standard momentum, which for photons is p u = (1,0,0, 1). Accordingly, standard 
right and left circular polarization vectors are = (1, ±i, 0)/y / 2. If we are interested in 
linear polarization, all we have to do is to use Eq. (13) with a + = («_)*, so that a(k) is 
real. In general, a(k) corresponds to elliptic polarization. 

Under a Lorentz transformation A, these states become \h^,a(k^)), where kA is the 
spatial part of a four-vector k\ = Ak, and the new polarization vector can be obtained by 
an appropriate rotation that is described below. For each k a polarization basis is labeled 
by the helicity vectors, 



e 



k 



i2(k)6± (15) 



and the corresponding quantum states are just |k, e k ). As usual, k denotes the unit 3- 
vector in the direction of k. The matrix that rotates the standard direction (0,0, 1) to 
k = (sin 9 cos 0, sin 9 sin 0, cos 9) is 



R(k) = 



( cos 9 cos — sin cos sin 9 \ 

cos 9 sin cos sin sin 9 
\ -sin 6* cos 9 J 



(16) 



and likewise for k A . 

Under a general Lorentz transformation, be it a rotation or a boost, helicity is preserved, 
but the states (and corresponding geometric vectors) acquire helicity-dependent phases, 

a + e+ + «_e k - a + e* A *e+ A + a.e'^e^, (17) 

where the explicit expressions for £(A, k) are given in [6] The superselection rule that was 
mentioned above makes it impossible to define a reduced density matrix in the usual way. 
We can however define an "effective" reduced density matrix for polarization, as follows. 
The labelling of polarization states by Euclidean vectors e£, and the fact that photons 
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are spin-1 particles, suggest the use of a 3 x 3 matrix with entries labelled x, y and z. 
Classically, they correspond to different directions of the electric field. For example, when 
k = k A z, only p xx , p xy , p yy are non-zero. 

For a generic photon state let us try to construct a reduced density matrix p xx that 
gives the expectation value of an operator representing the polarization in the x direction, 
irrespective of the particle's momentum. To have a momentum-independent polarization 
is to tacitly admit longitudinal photons. Therefore, in terms of real transversal photons 
of a momentum k, this problem suggests to find the direction that is perpendicular to k 
and closest to an arbitrary unit vector x. That is, we are looking for a unit Euclidean 
complex vector e x (k) such that k • e x (k) =0 and x • e x (k) = max. 

Momentum-independent polarization states thus consist of physical (transversal) and 
unphysical (longitudinal) parts with a polarization vector e l = k. For example, a polar- 
ization state along the x-axis is 

|x>=x + (k)|ej>+x_(k)|ei:>+x,(k)|4>, (18) 

where x±(k) = e k -x, and a^(k) = xk = sin#cos0. It follows that |x + | 2 +|x_| 2 +|:Ee| 2 = 1, 
and 

e x (k) = " +(k) ; k+ + "- (k)6k . (19) 

sjx\ + X^L 

Note that (x|y) = x • y = 0, whence 

|x)(x| + |y)(y| + |z)(z| = l (20) 
To the direction x corresponds a projection operator 

P xx = |x)(x| <g> l p = |x)(x| <g> J d/i(k)|k)(k|, (21) 

where l p is the unit operator in momentum space. The action of P xx on follows from 
Eq. (18) and (ekl e k) = 0- Only the transversal part of |x) appears in the expectation 
value: 

(*|Px*|*> = / rf/i(k)|/(k)| 2 |x + (k)«;(k)+x_(k)«l(k)| 2 . (22) 
Define the transversal part of |x): 

|b x (k)> = (|e+)(e+| + |e k )(e k |)|x) = s + (k)|e£> + x_(k)|e k ). (23) 

Likewise define |b y (k)) and |b z (k)). These three state vectors are neither of unit length 
nor mutually orthogonal. For k = (sin 9 cos 4>, sin 9 sin 0, cos 9) we have 

|b x (k)) = [(cos # cos + ?sin0)|e k ) + (cos # cos — ?sin0)|e k )]/v / 2, (24) 
= c(0,0)|k,e x (k)), (25) 



where e x (k) is given by Eq. (19), and c(9, <fi) = Jx 2 + + x 2 
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Finally, a POVM element E xx which is the physical part of P xx , namely is equivalent 
to P xx for physical states (without longitudinal photons) is 

E xx = J dfi(k)\k,b x (k))(Kb x (k)\, (26) 

and likewise for other directions. The operators E xx , E yy and E zz indeed form a POVM in 
the space of physical states, owing to Eq. (20). It then follows from Eq. (23) and similar 
definitions for other directions that, for any k, 

M^XMk)! + |b w (k))(b w (k)| + Mk))(b,(k)| = i ±k , (27) 

where is the identity operator in the subspace of polarizations orthogonal to k. 

To complete the construction of the density matrix, we introduce additional directions. 
Following a standard practice of state reconstruction [7], we consider P x+Z ^ x+Z , P x ~i Z ,x-iz 
and similar combinations. For example, 

P x+ z,x+z = |(|x) + |z»((x| + (z|) <g> t p (28) 
The diagonal elements of the new polarization density matrix are defined as 

p mm = (&\E mm \V) , m = x,y,z, (29) 
and the off-diagonal elements are recovered by combinations such as 

p xz = (tf(|x)(z| <g> = (V\E X+Z}X+Z + E x _ iZtX _ iz - E xx - E ZZ \V), (30) 

where we denote |x)(z| ® t p as P xz , and by E xz its "physical" part. We then get a simple 
expression for the reduced density matrix corresponding to the polarization state |a(k)): 

Pmn = (*|£m„|*) = / dfjlQz) \ f (k) | 2 (tX (k) | b m (k) ) (b n (k) | OL (k) ) , TO, 71, = X, y, Z. (31) 

Our basis states |k, e^) are direct products of momentum and polarization. Owing to 
the transversality requirement e(k) • k = 0, they remain direct products under Lorentz 
transformations. All the other states have their polarization and momentum degrees 
of freedom entangled. As a result, if one is restricted to polarization measurements as 
described by the above POVM, there do not exist two orthogonal polarization states. In 
general, any measurement procedure with finite momentum sensitivity will lead to the 
errors in identification. This can be seen as follows. 

Let two states |<&) and |^) be of the form in Eq. (12). Their reduced polarization den- 
sity matrices, p$ and p^, respectively, are calculated using Eq. (31). Since the states are 
entangled, the von Neumann entropies of the reduced density matrices, S = — tr(plnp), 
are positive [8]. Therefore, both matrices are at least of rank two. Since the overall 
dimension is 3, it follows that tr(p$p*) > and these states are not perfectly distinguish- 
able. An immediate corollary is that photon polarization states cannot be cloned perfectly. 
This is because no-cloning theorem, in its various versions [9], forbids an exact copying 
of unknown non-orthogonal states. 

To quantify the distinguishability of a pair of quantum states, we shall use the simplest 
criterion, namely the probability of error Pe, defined as follows: an observer receives a 
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single copy of one of the two known states and performs any operation permitted by 
quantum theory, in order to decide which state was supplied. The probability of a wrong 
answer for an optimal measurement is [10] 

P E (pi,P2) = h + itT\p 1 -p 2 \, (32) 



where, for any operator O, the operator \0\ is defined as y(0^0). As shown below, the 
distinguishability of polarization density matrices depends on the observers' motion. First 
we present some general considerations and then illustrate them with a simple example. 

Let us take the z-axis to coincide with the average direction of propagation so that 
the mean photon momentum is fc^z. Typically, the spread in momentum is small, but 
not necessarily equal in all directions. Usually the intensity profile of laser beams has 
cylindrical symmetry, and we may assume that A x ~ A y ~ A r where the index r means 
radial. We may also assume that A r ^> A z , since the momentum spread along the 
average propagation direction does not necessarily increase the mixness of the reduced 
density matrix, as we shell see below. We then have 

/(k) oc h[{k z - k A )/A z ] f 2 (k r /A r ). (33) 

We approximate 

6 ta,n6 = k r /k z ps k r /k A . (34) 

In pictorial language, polarization planes for different momenta are tilted by angles up 
to ~ A r /kA, so that we expect an error probability of the order A 2 r jk\. In the density 
matrix p mn all the elements of the form p mz should vanish when A r — > 0. Moreover, 
if A z — > 0, the non-vanishing xy block goes to the usual (monochromatic) polarization 
density matrix, 

/ M 2 (3 \ 

Ppurc =\ (3* 1 - \a\ 2 . (35) 
\ / 

As an example we consider two states which, if the momentum spread could be ignored, 
would be \k A z : e k A z)- To simplify the calculations we assume a Gaussian distribution: 

f{ y )=Ne -^k A? /2^ e -ki/2A^ (36) 

where iV is a normalization factor and A z <C A r . Moreover, we take the polarization 
components to be = R(k)e^ . That means we have to analyze the states 

|*±> = |^(k)/(k)|e± k), (37) 

where /(k) is given above. 

We expand -R(k) up to second order in 6. Reduced density matrices are calculated by 
techniques similar to those for massive particles [3], using rotational symmetry around 
the z-axis and normalization requirements. The leading order in Q = A r /k A gives 

(38) 
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and p_ = p* + . At the same level of precision, 

P E (p + ,p.) = A 2 r /4k 2 A . (39) 

It is interesting to note that an optimal strategy for distinguishing between these two 
states is a polarization measurement in the xy-plane. Then the effective 2x2 density 
matrices are perfectly distinguishable, but there is a probability Q 2 /2 that no photon be 
detected at all. The above result was valid due to the special form of the states that we 
chose. Potential errors in the upper 2x2 blocks were averaged out in the integration 
over 0. 

We now calculate Bob's reduced density matrix. We again assume that Bob moves 
along the z-axis with a velocity v. Recall that reduced density matrices have no transfor- 
mation law [3], (only the complete density matrix has one) except in the limiting case of 
sharp momenta. The only way to calculate Bob's reduced density matrix is to transform 
the complete state, and only then take a partial trace. According to Eqs. (31) and (37), 
reduced density matrices in both frames are given by the expression 

(p±)mn = J rf/i(k)|/(k)| 2 (i?(k)e p ± |b m (k))(b n (k)| J R(k)e p ± ). (40) 

Note that pure boosts preserve the orientation of the coordinate axes in 3-space, and 
therefore do not affect the indices of p mn . The measure /i(k) is Lorentz-invariant and 
€p are constant by definition. Since / is a scalar function, it transforms as /'(k) = 
/(k A -i), where primes indicate Bob's frame, as in Eqs. (6-9). This is the only frame 
dependent expression in (40), since the phase factors e ±l ^ cancel out for the helicity 
eigenstates. Therefore, there are two equivalent methods to calculate Bob's polarization 
density matrix. One is to change the argument of /(k) to k A i, and another is to change 
the argument of the rotation matrix -R(k) to k A . Using the second option, we obtain 

(p' ± ) mn = j rf/i(k)|/(k)| 2 ( J R(k A )e p ± |b m (k))(b n (k)| J R(k A )e p t ). (41) 

A boost along the z-axis preserves k r and <fi. On the other hand, from Eq. (7) it follows 
that 



^Wirf (42) 

Thus, at leading order in 9 we have 



e>*J±±^k r /k A *J^e, (43) 

V 1 — v f 1-u 

which is substituted into i?(k A ). Since everything else in the integral remains the same, 
the effect of relative motion is given by a substitution 



n-^n. (44) 

It follows that 

P' E = \^Pe, (45) 



8 



which may be either larger or smaller than Pe- As expected, we obtain for one-photon 
states the same Doppler effect as in the preceding classical calculations. 

Although reduced polarization density matrices have no general transformation rule, 
the above results, as well as the analysis of massive particles [3], show that such rules 
can be derived for particular classes of experimental procedures. We can then ask how 
these effective transformation rules p' = T[p] fit into the framework of general state 
transformations. A general state transformation T is usually required to be completely 
positive (CP), namely [11, 12], 

T[p]=Y,M lP Ml (46) 

i 

where the Mj are bounded operators. It can be proved that distinguishability, as expressed 
by natural measures like cannot be improved by any CP transformation [10]. It is also 
known that the CP requirement may fail if there is a prior entanglement of p with another 
system [13]. Since from [3] and Eq. (45) it follows that in our case distinguishability can 
be improved, we conclude that these transformations are not completely positive. The 
reason is that the Lorentz transformation acts not only on the "interesting" polarization 
variables, but also on the "hidden" momentum variables that we elected to ignore and to 
trace out. 

This technicality has one important consequence. In quantum information theory 
quantum channels are described as completely positive maps [14, 15, 16] that act on 
qubit states. Qubits themselves are realized as particles' discrete degrees of freedom. 
If relativistic motion is important, then not only does the vacuum behave as a noisy 
quantum channel, but the very representation of a channel by a CP map fails. 
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